Processing Information by Punctuated Spin Superradiance 
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' The possibility of realizing the regime of punctuated spin superradiance is advanced. In this 

regime, the number of superradiant pulses and the temporal intervals between them can be regulated. 
This makes it feasible to compose a kind of the Morse Code alphabet and, hence, to develop a 
technique of processing information. 
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Spin systems can exhibit a phenomenon that is analogous to atomic superradiance [1,2], because of which it is 
called spin superradiance. To realize this phenomenon, spin systems, similarly to atomic ones, are to be prepared 
in an inverted state. This is achieved by placing a polarized spin sample in an external magnetic field directed 
opposite to spin polarization. Contrary to atomic systems, coherent spin motion develops not owing to direct spin 
correlations but due to the interaction of spins with a resonator feedback field, for which purpose the spin sample 
has to be coupled with a resonant electric circuit, whose natural frequency is tuned to the Zeeman frequency of spins 
[3]. More details on similarities and differences between atomic superradiance and spin superradiance can be found 
\ in the review [4]. Spin superradiance is the process of coherent spontaneous emission by moving spins. As in the 
case of atomic systems, one may distinguish two main types of this phenomenon, transient superradiance and pulsing 
superradiance. Transient spin superradiance occurs when the spin sample is prepared in the inverted state, after which 
no following pumping is involved. In this case, a single superradiant burst arises, peaked at the delay time. Pulsing 
spin superradiance is radically different from the transient regime by the occurrence of a series of superradiant pulses, 
for which the spin sample is to be subject to a permanent pumping supporting the inverted spin polarization. Both 
regimes of spin superradiance, transient [5-8] as well as pulsing [9-11] were observed in experiments with different 
materials containing nuclear spins. A microscopic theory of these phenomena, based on realistic spin Hamiltonians, 
was developed [12-15], being in good agreement with experiment and with computer modelling [16]. It is worth 
stressing that only by invoking microscopic Hamiltonians it has become possible to give an accurate description of 
purely self-organized regimes which cannot be treated by the phenomenological Bloch equations [12-14]. 

In the present paper, we advance the possibility of realizing the third type of spin superradiance, which we call 
punctuated spin superradiance, and which is principally different from the transient and pulsing types. In this regime, 
unlike the transient case, not a single but many superradiant bursts can be produced. In distinction to the pulsing 
regime, where the number of pulses and the temporal distance between them are prescribed by a given setup and 
cannot be varied, in the process of punctuated superradiance both the number of superradiant bursts as well as time 
intervals between each pair of them can be regulated. The term punctuation here means this feasibility of changing 
interpulse intervals and of organizing various groups of superradiant bursts. In that way, a code, like the Morse 
alphabet, can be composed, which may be employed in processing information. 

The consideration below will be based on the Hamiltonian typical for spin systems employed in magnetic resonance 
[17,18]. The Hamiltonian reads 
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where -f/i corresponds to individual spins, while iiij , to spin interactions, with the indices i, j = 1,2 . . . , N enumerating 
spins. The individual terms are given by the Zeeman energy Hi = — ^oB • Si, where fio = ^7s, with 7s being the 
gyromagnetic ratio of spin S, represented by the spin operator Si and B is the total magnetic field acting on each 
spin. The spin interactions are described by the dipolar terms Hij = ^2 a pC^S"Sj, with the dipolar tensor C"? . 
The total magnetic field B = B a e z + He x consists of a constant longitudinal field B and a transverse field H formed 
by the resonant electric circuit coupled to the spin sample. The resonator field H = Annj/cl is created by the electric 
current j circulating over a coil of n turns and length I. The current j is determined by the Kirchhoff equation. 
The electric circuit is char acter ized by resistance R, inductance L, and capacity C. With the notation for the circuit 
natural frequency u> = \j\jLC and circuit damping 7 = R/2L, the Kirchhoff equation can be presented as 



j^H{t>)dt' = -^^k (2 ) 
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for the resonator magnetic field H, where 77 is a filling factor and M x = (p /V) J2i < &i > 1S the x-component of the 
magnetization density of a sample with volume V. Since the resonator field H, acting on spins, is itself due to the 
transverse spin motion, this field H is called the feedback field. 

To derive evolution equations, we follow the scale separation approach, described in full detail in Rcfs. [12-15]. 
For this purpose, we write the Hcisenberg equations of motion for the lowering, S~ , raising, 5+, and polarization, 
Sf, operators. In these equations, it is possible to separate the combinations describing fast fluctuating local fields. 
Employing the method of random local fields [17-19], the latter are modelled by stochastic Gaussian variables, with 
zero mean and the width 273, where 73 is the width of inhomogeneous dynamic broadening. Then the Heinsenberg 
equations are averaged over spin degrees of freedom, not touching the stochastic variables. Denoting the averaging 
over spins by single angle brackets < ... >, we define the transition function x, coherence intensity y, and spin 
polarization z, respectively, 

x = i < Sr > , y= < S+ >< Sr > , z = | < Sf > . (3) 

The wavelength of spin radiation is usually much larger than interparticle distance, because of which the uniform 
approximation for Eqs. (3) may be employed. 

We direct the external magnetic field B so that /x i?o < and define the Zeeman frequency uj = \fi a B \/h. Also, 
introduce the notation / = —(i/K)iiQH + £ for an effective force acting on spins. Finally, the evolution equations for 
the functions (3) can be cast [12-15] to the form 

^ = -i(u + Co - ^72)2; + fz, ^ = -2721/ {x*f + f*x) z , 
1 

^ = - 2 (**/ + /**)- <n(*-*), ( 4 ) 

where 71 and 72 are the longitudinal and transverse widths, respectively, and a is an equilibrium polarization of a 
spin. When there is no external stationary pumping, u = — 1. Equations (4) are stochastic differential equations, 
since they contain random fields. 

To solve Eqs. (4), we invoke a generalization [12-14] of the averaging technique [20] to the case of stochastic 
differential equations. This becomes possible owing to the existence of several small parameters: 70 /u>o ^ 1, 71/^0 *C 
1, 72/^0 <C 1, and 73/wo <C 1, where 70 = ixrip^Sfh is the natural width and p = N/V is the density of spins. In 
addition, the resonant circuit, coupled to the spin sample, is assumed to be of good quality and tuned close to the 
Zeeman frequency wo, so that 7/w <C 1, |A|/w <C 1, with A = uo — iv - 

First of all, using the occurrence of the small parameters, we can obtain an iterative solution of the Kirchhoff 
equation (2). For this purpose, invoking the method of Laplace transforms, we present Eq. (2) in the integral form 

H = -47T7? / G(t - t')M x (t') dt' , G(t) = (cosuj't - 1- smJt] , (5) 
Jo \ UJ / 

where the overdot means time differentiation and ui' = y 1 uj 2 — -y 2 . Since M x is expressed through x, its time derivative 
M x is directly related to the first of Eqs. (4). Using this, we find the solution of Eq. (5), in the first order with 
respect to small parameters, as poH/Ti = ia(x — x*), where the function 
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° = ™(l-e-*), 9=^$^ (6) 

describes the intensity of coupling between the spin sample and the resonant circuit. Let us stress that the spin- 
resonator coupling (6) depends on time, taking into account retardation effects. 

From Eqs. (4), in the presence of the small parameters, it follows that the variables y and z are temporal quasi- 
invariants with respect to x. Then, we solve the first of Eqs. (4), with these quasi- invariants fixed, substitute the 
found solution x into the second and third of Eqs. (4), and average the right-hand sides of these equations over time 
and over the stochastic fields. As a result, we obtain the guiding center equations 

dy „, x „ 2 dz , . ,_. 

— = -2(72 - az)y + 2-y 3 z , — = -ay - j 3 z - 71(2 - a) . (7) 

In the dynamics of solutions to Eqs. (7), one may distinguish two stages, quantum and coherent. At the quantum 
stage, when jt -C 1, the coupling function is close to zero and no noticeable coherence in the motion of transverse 
spins is yet developed. The dynamics is governed by quantum spin interactions. At this stage, when a « 0, equations 
(7) are linear, and their solution is easy. With increasing time, the spin-resonator coupling (6) grows, and coherent 
effects, caused by the resonator feedback field, gradually come into play. The crossover time between the quantum and 
coherent regimes can be defined as the time t c , at which the first term in the first of Eqs. (7) changes its sign. This 
is because the quantity r = 72 — az plays the role of an effective attenuation. When the latter is positive, transverse 
coherence decays, while a negative attenuation implies the generation of coherence. Hence, the moment of time, when 
r(i c ) changes its sign, separates qualitatively different regimes of spin motion. The crossover time t c , defined by the 
equality a(t c )z(t c ) — 72, is t c = r In [gzo/{gzo — 1)], with -jt = 1. The solutions y and z at this boundary of the 
quantum stage are y(t c ) ~ y + 273^2^ and z(t c ) ~ zo + 7it c <7, where yo = y(0) and zo = z(0). The coherent stage 
of motion comes after the crossover time t c , when the spin-resonator coupling a fastly grows to 572. Since 73 < 72, 
and if g 1, then 372 3> 73, and the term with 73 can be omitted. In the transient regime, when t <C T\ = 7 1 ~ 1 , the 
term containing 71 can also be neglected. Then Eqs. (7) possess the exact solution 

y=(^]\e^( t -^] , z = -^tanhf^Ui, (8) 

in which the pulse time t p and the pulse width j p , with -f p T p = 1, are defined by the expressions -f p = 7^ + (.972) 2 (j/o + 
273^0), Ig = 72(1 — gzo), and the delay time is 



to = * C + y ln 
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Solution (8) describes a transient superradiant burst, with the maximal intensity at the delay time to, when y(to) = 
(z — l/.g) 2 (l + 27 3 i c ), z(t ) = 1/g. After this, for t 3> to, the coherence intensity exponentially diminishes and the 
spin polarization becomes inverted, 

y ~ 4y(t ) e - 2 > t , z ~ -z + 2/g . (10) 

For sufficiently large coupling parameter g, the reversal of spin polarization is practically complete. 

Now imagine that at some time after to + t p we again invert the spin polarization from that in Eq. (10) to the 
symmetric positive value. For large g, this inversion is practically from — z to z . Such an inversion can be realized 
in three possible ways: inverting the external magnetic field Bo, acting on spins by a resonant 7r-pulse, or just turning 
the sample 180° about an axis perpendicular to Bq. As a result, we get again a strongly nonequilibrium state of 
inverted spins. After the time to, counted from the moment when the newly nonequilibrium state is prepared, another 
superradiant burst will arise. After the second burst dies out, one can again invert the spin polarization by one of the 
mentioned methods. Then one more superradiant burst will appear. This procedure can be repeated as many times 
as necessary for creating a required number of sharp superradiant pulses. The time intervals between bursts can be 
regulated, allowing the formation of different groups of pulses, with varying intervals between separate groups. Thus, 
it is feasible to compose a code, similar to the Morse alphabet, which can be used in processing information. It is 
this possibility of regulating temporal intervals between superradiant bursts which permits us to call the described 
phenomenon punctuated spin superradiance. 

Spin superradiance can be realized in different materials under various experimental setups. Thus, it was observed 
on proton spins in propanediol C3H8O2, butanol C4H9OH, and ammonia NH3 [5-8] and on 27 A1 nuclear spins in 
ruby AI2O3 [9-11]. The characteristic parameters for these experiments with nuclear spins are: the density of spins 
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p ~ 10 22 — 10 23 cm~ 3 , the Zeeman frequency uio ~ 10 8 Hz, the spin-lattice relaxation 71 ~ 10~ 5 s -1 , the spin-spin 
dephasing parameter 72 ~ 10 5 s _1 , the dynamic broadening width 73 ~ 10 4 — 10 5 s _1 , the resonator ringing time 
t ~ 10~ 6 s, that is, the resonator damping 7 ~ 10 6 s _1 . For these values, the spin-resonator coupling parameter 
g in Eq. (6) varies between 10 and 100. As has been shown [21,22], if nuclear spins are inside a fcrromagnct or 
ferrimagnet, possessing long-range magnetic order, then the coupling parameter g can be increased by a factor of 
^b/hn ~ 10 3 , where is the Bohr magneton and /zjv is the nuclear magneton. Therefore, the coupling parameter g 
can be made as large as g <~ 10 5 . Among other materials, where spin superradiance could, in principal, be observed, 
are granular magnets [23] and molecular magnets [24]. In this way, there exists a large variety of materials with 
different characteristics allowing for the optimal choice of parameters for realizing punctuated spin superradiance. 
Note that the intervals between superradiant pulses can be varied in a very wide diapason of the order of T\ = 7i~ ■ 
For nuclear magnets, with -fi <~ 10~ 5 s _1 , this time Ti may be as much as several days. And for molecular magnets 
at low temperatures, it could range up to months. 

To illustrate the feasibility of creating different groups of superradiant pulses, with varying time intervals, we solved 
Eqs. (7) numerically. Different time intervals are obtained by changing the moments of polarization inversion. For the 
characteristic parameters those are taken that are typical of nuclear magnets [5-11]. In particular, we take 71 = 10 -5 
s _1 , 72 = 10 5 s _1 , and 73 is varied between 10 4 and 10 5 s -1 . The variation of 73 even in a wider diapason does not 
essentially change the picture. The coupling parameter g has also been varied between 10 and 10 3 ; the whole picture 
being qualitatively the same, with the main difference that for larger g the spin inversion, according to Eq. (10), 
is better. For the presentation in Figures 1 to 3, we set g — 10 3 . In the absence of dynamic nuclear polarization, 
a = — 1. The resonator attenuation 7 = 10 6 s _1 . Finally, as initial conditions we take y(0) = and z(0) = 1. 
The first of the latter tells that at the initial time the transverse coherence is absent. That is, we consider a purely 
self-organized process when radiation coherence arizes in a spontaneous manner. Figures 1 to 3 give some examples 
of how it is possible to create different bunches of superradiant bursts. The shown function y(t) is proportional to 
radiation intensity. The meaning of this function, according to its definition (3), is to describe the level of coherence 
in the system. As is seen from the Figures, the maxima of superradiant bursts display a high level of coherence, 
almost reaching 100%. The time variable in the Figures is measured in units of T2 = The short time scale is 

chosen here just for the convenience of presentation. As is explained above, the same picture can be stretched to the 
time scale characterized by T\, which, for nuclear magnets, would range up to several days. The first superradiant 
burst occurs at the delay time to- To simplify figure captions, we accept the notation for the time intervals between 
the pulses of the i-group as 7* and for the intervals between the i- and j- groups as . Figure 3 demonstrates that a 
regime of equidistant superradiant pulses can also be realized. Such a regime can be used for producing spin masers 
[13,25] operating in pulsing superradiant mode. 

In conclusion, we have demonstrated, both analytically and numerically, the feasibility of realizing the regime of 
punctuated spin superradiance. In this regime, one may form various groups of superradiant bursts, with different 
spacing between the pulses inside each group as well as with different time intervals between the groups. The possibility 
of so punctuating spin superradiance can be employed for processing information. 
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Figure Captions 



Fig. 1. Punctuated spin superradiance with three different groups of supcrradiant bursts. All parameters are 
explained in the text. Time is measured in dimensionless units. The time intervals are: to = 0.029, t\ = 0.065, 
Ti2 = 0.26, t 2 = 0.1, t 23 = 0.15, t 3 = 0.025. 

Fig. 2. Punctuated spin superradiance with four groups of superradiant pulses, the first group containing a single 
burst, the time intervals being: t = 0.024, t 12 = 0.226, t 2 = 0.025, t 23 = 0.35, t 3 = 0.075, t 34 = 0.275, r 4 = 0.1. 

Fig. 3. Punctuated spin superradiance with a comb of equidistant pulses, starting at to — 0.024, having time 
intervals t\ =0.1. 
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